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Abstract — This paper is concerned with the error performance 
analysis of binary differential phase shift keying with differential 
detection over the nonselective, Rayleigh fading channel with com- 
bining diversity reception. Space antenna diversity reception is as- 
sumed. The diversity branches are independent, but have noniden- 
tically distributed statistics. The fading process in each branch is 
assumed to have an arbitrary Doppler spectrum with arbitrary 
Doppler bandwidth. Both optimum diversity reception and sub- 
optimum diversity reception are considered. Results available pre- 
viously apply only to the case of first and second-order diversity. 
Our results are more general in that the order of diversity is arbi- 
trary. Moreover, the bit error probability (BEP) result is obtained 
in an exact, closed-form expression which shows the behavior of 
the BEP as an explict function of the one-bit-interval fading corre- 
lation coefficient at the matched filter output, the mean signal-to- 
noise ratio per bit per branch and the order of diversity. A simple, 
more easily computable Chernoff bound to the BEP of the opti- 
mum diversity detector is also derived. 



I. INTRODUCTION 

An efficient modulation technique for communication over a 
nonselective Rayleigh fading channel is differential phase shift 
keying (DPSK) with differential detection and combining di- 
versity reception at the receiver [1]. Most previous analyses 
assume that the fading processes on the diversity branches are 
both independent and identically distributed (i.i.d.) [1] — [9]. 
However, in practical systems, this assumption is not always 
true [10]. The mean-square values of the fading attenuations 
are usually not identical for all diversity branches [11], [12], 
i.e., the average signal-to-noise ratio (SNR) of each branch is 
distinct. Therefore, it is of great practical interest and impor- 
tance to study the performance of combining diversity reception 
for differentially detected DPSK over independent, nonidenti- 
cal Rayleigh fading channels. With nonidentical fading chan- 
nels, the combiner described in [7] and [8] is suboptimal, and 
it is well known that the BEP performance can be improved 
if each branch differential detector output is weighted, before 
combination, according to its channel condition [13]. We re- 
fer to this combining method with optimum weights for the 
branches as optimum combining reception. The bit error proba- 
bility (BEP) performance of optimum combining reception was 
analyzed in [14]. However, the results are limited to the case 
of second-order diversity. Moreover, to the best of the authors' 
knowledge, as far as the exact, explict, closed-form BEP ex- 
pression is concerned, there has not been much progress. The 
main contribution of this paper lies in providing the exact anal- 



ysis of the BEP of binary DPSK with arbitrary order of diver- 
sity, both for the optimum combining reception and the subopti- 
mum combining reception. The fading process is assumed here 
to have an arbitrary Doppler spectrum with arbitrary Doppler 
bandwidth. The BEP results here are derived in exact, closed- 
form expressions which depend explictly on the SNR per bit per 
branch, the fading correlation coefficient at the matched filter 
output over a bit interval, and the order of diversity. The solu- 
tions do not require any numerical integration for their actual 
evaluation. Moreover, a simple, tight, more easily computable 
Chernoff upper bound for the optimum diversity detector case 
is also obtained. 

The paper is organized as follows. In Section II, the system 
model and assumptions are introduced. The analyses for the 
BEP and the upper bound are given in Section III. Section IV 
presents numerical results. 

II. SYSTEM MODEL AND ASSUMPTIONS 

With space antenna diversity reception over L independent, 
nonidentical, frequency nonselective, Rayleigh fading branches 
with additive, white Gaussian noise (AWGN), the received sig- 
nal over the ith branch, i = 1,2, • • ■ , L, during the fcth bit 
interval [kT, (fc + 1)T) is given, after matched filtering and 
sampling at time t = (k + 1)T, by the statistic Sj(fc), where 



z i (k)=E 1 b /2 e^~a t (k) + u l (k). 



(1) 



The overhead ~ denotes a complex quantity, and the superscript 
* denotes the complex conjugate. Eb is the energy per bit, T is 
the bit duration, 0(fc) is the data-modulated phase for the fcth 
bit interval. The phase transition, A<fi(k) = <fi(k) — <p(k — 
1), between the fcth and the (fc — l)th bit intervals carries the 
signal information. For binary DPSK, the data 0-bit is mapped 
onto the phase change A<p(k) = 0, and the data 1-bit onto the 
phase change A<j>(k) = it. We assume that the two data bits are 
equally likely. The multiplicative distortion fi;(fc) is given by 



Oi(fc) 



(fc+l)T - 



kT 



Cj(t) 

T 



dt. 



(2) 



Here, {cj(t)}-j is a set of independent, nonidentically dis- 
tributed, lowpass, complex Gaussian random processes with 
E [5i(t)] = and E[5;(t)5*(t - t)} = 2i^(r) for each i. 
Each Cj (t) represents the complex gain due to frequency non- 
selective Rayleigh fading of the ith branch. We assume that 



the spectrum of each Ci(t) is even so that the in-phase compo- 
nent Re [ci(t)} and the quadrature phase component Im [<k(t)] 
are independent with the same co variance function J?j(r). In 
(1), a rectangular data pulse shape g(t), where g{t) — 1/y/T 
for < t < T and zero elsewhere, is assumed so that each 
matched filter has a rectangular low-pass-equivalent impulse re- 
sponse hi(t) = g(T — t) for all i. This leads to the expression 
for 5j(fc) in (2). Thus, {5,i(k)} k is a sequence of zero-mean 
complex Gaussian random variables with covariance function 



A. BEP analysis 

First, we consider the case of optimum reception. Because 
the BEP is the same whether a 0-bit or a 1-bit is sent, we assume 
that A0(fc) = 0. Without loss of generality, it is assumed that 
4>{k) = 4>(k - 1) = 0, and the BEP P b is given by 



P b = P < Rc 



Rid) 



i 



E[oi{k)^(k-j)] 



1 

2^2 



{k+l)T r (k-] + l)T 



kT 



(k-j)T 



Ri(t - T)dTdt. (3) 



The filtered noise u%{k) is given by 

f (k+1)T n t (t) 
IkT Vt 



i>i{k) 



dt. 



(4) 



Here, {fii(t)} i=1 is a set of i.i.d., lowpass, complex AWGN 
processes with E [n,i(t)] = and E[fii(t)h*{t - t)] = N S(t) 
so that {vi(k)} k is a sequence of i.i.d., zero-mean, complex 
Gaussian variables with covariance function for each i 



E[P i (k)P*(j))=N 6 kj 



(5) 



where S k j is Kronecker delta function. 

For each i , Ci (t) and n, (t) are mutually independent. For 
i 3i {ci(t) , fii(t)} are independent of {cj(t),hj(t}}. The 
channel branches are nonidentical since the covariance func- 
tion Ri(r) depends on i. For convenience of later application, 
the following parameters are defined. The quantity pi = ^' 
is the fading correlation coefficient at the matched filter out- 
put over an interval of T, and it is a measure of the fluctuation 
rate of the channel fading process; 7; = 2Eb ^^°> [ s the mean 
received SNR per bit over the ith branch. 

III. THE PERFORMANCE ANALYSIS 

In this section, we will derive exact, closed-form BEP ex- 
pressions and Chernoff bounds for differentially detected bi- 
nary DPSK for the optimum combining diversity detector and 
the suboptimum combining diversity detector. In the Appendix, 
we show that the optimum combining detector with differential 
detection makes its decision on the fcth transmitted bit using the 
likelihood ratio test 



Re 



^ w l Si(k) z*(k - 1) 



(=1 



bit 

I 

1 bit 



where 



pni 



(6) 



(7) 



(l+ 7i ) 2 -(ft70 2 ' 
In the case of suboptimum combining diversity detection, one 
has Wi = 1, for i = 1, 2, • • • , L. Then the likelihood ratio test 
becomes 



Re 



i=l 



bit 

% 

1 bit 



(8) 



This is the same as the optimum combining diversity detection 
in the i.i.d. channel case. 



Wi Zi(k) z*(k - 1) 



i=l 



< 



A<j>(k) = 



(9) 



P 



Wi\zi{k) + Zi(k -i)i < Wi\zi(k) ~ Zi(k - 1)| 



The alternative case where the 0-bit sent is due to 4>{k) = <fi(k— 
1) = 7r gives an identical result. With 4>{k) = <fi(k — 1) = 
0, one can see from (1) that for each i = 1, 2, • • • , L, Zi{k) 
and Zi(k — 1) are both zero-mean complex Gaussian random 
variables (since each one is a sum of two zero-mean complex 
Gaussian random variables) with variances 

E[\Uk)\ 2 ] = E[\zi{k - 1)| 2 ] = 2^(0) + N . (10) 

Therefore, the quantity Zj(fc) + Zi(k — 1) is complex Gaussian 
with mean zero and variance 4Ef,i?j(0) + AE b Ri{l) + 2N , 
and 5i(k) — Zi(k — 1) is complex Gaussian with mean zero and 
variance AEbRi (0)—4EbRi ( 1 ) + 2N . Moreover, the quantities 
Zi(k) + Zi(k — 1) and 2,(fc) — Zi(k — 1) are uncorrected since 

E{[ Zl (k) + Zi(k - 1)] [Zi(k) - Zi (k - 1)}*} = E[\zi(k)\ 2 ] 
-E[\z i (k-l)\ 2 ]-j2E{Im[z i {k)z*(k-l))} = (11) 

where j — s/—l. The last step follows from (10) and the result 
E[zi(k)z*(k — 1)] = 2EbRi(l), which is a purely real quantity. 
Thus, Zi(k) + Zi(k — 1) and Zi(k) — Zi(k — 1) are independent. 
To proceed with computing (9), define 



Xi = Wi\z.i(k) + li(k 
and 



yi = Wi\zi(k) - Zi{k - l)Y 



i=l i=l 

Then the BEP Pb in (9) can be rewritten as 



P b = P ( X < Y ) 



(12) 
(13) 

(14) 



We note that each Xi and each in (12) are chi-square- 
distributed with two degrees of freedom. The probability den- 
sity function (PDF) can be evaluated as [1, (2.1-126)] 



PiVi) 



1 



Wi[4E s Ri{0) + AE s R l {\) + 2N ] 



(15) 



exp 



^[4^,^(0) + 4E s Ri(l) + 2N ] 
1 



Wi[4E a Ri(0) -AE a Ri{l) 
Vi 



2Nn 



(16) 



exp 



w i [4E s R i (0)-4E s R i (l) + 2N 



For each i , x% and yi are mutually independent, because zt (k) + 
Zi(k — 1) is independent of Zi(k) — Zi(k — 1). Moreover, due 



to the independent channel assumption, for i ^ j, {xt, yi} are 
independent of {xj,yj}. Therefore, the quantity X in (13) is 
independent of Y, The PDFs of X and Y can be obtained by 
using [1, (14.5-26)] and, after simplification, are given by 

L 



E 



A, 



lWi[4E a Ri(0)+4E a Ri{l) 
X 



2Nn 



x (17) 



p(Y) 



exp 



L 



Wi[4E B Rt{0) + AE a R t (l) + 2N Q ] 
Bi 



^ Wi[4E a Ri(p) - 4E a Ri(l) 
Y 



2N ] 



x (18) 



exp 



L Wi[AE a Ri(0) - AE a Ri(l) 
where Ai and Bi are defined as 



2A n 



A, 



with 



n 



B, = 



m—l,m^i 



n 



n— l,n^i 



Pi~Pn 



Pili 



Pi 



Pili 



1 



It + Pili 



1 + 7< - Pili 
Finally, the BEP in (14) can be evaluated as 

/>oo pY 

P b = P(X <Y)= / p(Y) / P (X)dX dY. 
Jo Jo 

Using (17) and (18) in (21), the BEP is obtained as 



L L 



Pi 



i=i j=i 



Pi 



(19) 



(20) 



(21) 



(22) 



Note that (22) gives an exact, explicit, closed-form BEP ex- 
pression. It depends only on the fading correlation coefficient 
Pi, the mean received SNR per bit per branch 7,, and the num- 
ber of diversity branches L. No numerical integration is needed 
for its actual evaluation. The quantities A, and Bi are given in 
(19), where on and Pi can be found in (20). 

Next, we consider the case of suboptimum reception. Again, 
it is assumed that <fi(k) = (f>(k — 1) = 0. The BEP of the subop- 
timum detector given in (8) can be obtained from the probability 



n 



p 



^2\zi(k) + Zi{k 

i=l 



1)| 2 <Y i \z i (k)-z i {k-l) 



i=l 



(23) 

It is apparent that when the weights w/s satisfy w; = 1 for 
all i = 1, 2, • • • , L, the probability expression (23) is identical 
to that given in (9). Consequently, the BEP performance of 
suboptimum combining diversity reception is identical to that 
given in (22), where Ai and Bi are defined in (19). It can be 
shown that for suboptimum reception oti and Pi are given by 



a, 



1+Ji + Pili, Pi = 1 + 7i _ Pili 



(24) 



Although exact and explicit, (22) is cumbersome and inconve- 
nient to use when L is large. Further, the behavior of the error 
probability (22) as a function of the various system parameters 
can not be easily seen. In the sequel, we will develop a simple, 
more easily computable Chernoff upper bound on the BEP. 



B. Chernoff bound analysis 

First, consider the optimum reception case. Applying the 
Chernoff bound [1, section 2.1.5] to (22) and noting that X is 
independent of Y, we have 



P h < E 



-s(X-Y) 



= E [er sX ] E [e sY ] 



(25) 



where s > is the parameter to be optimized. Using (12) and 
(13) in (25), and noting that for i ^ j, Xi is independent of Xj 
and yi is independent of yj, we have 



E [e- sX ] 



n4 



i=l 
L 



E [e-] = n4 



— sWi \Si{k)-\-Si(k— 1) I 



■ Wi\zi{k)-Zi{k~l)Y 



(26) 



(27) 



To proceed with computing (26) and (27), we use the following 
well-known lemma [15, (7.67)] 

Lemma 1: If x is a Gaussian random variable with mean v 
and variance a 2 and e is any complex constant with real part 
less than (2cr 2 ) -1 , then 



E 



1 



VI - 2ea 2 



;i/ 2 /(l-2eo- 2 ). 



Re(e) < A, 



Applying this lemma to (26) with e = —SWi, and (27) with 
e = SWi, an d noting that Si(k) + Zi(k — 1) and z,(fc) — Zi(k — 1) 
are zero-mean complex Gaussian with variances AEhRi(Q) + 
4E b R,(l)+2N and 4E b R l (0)-4E b R l (l)+2N , respectively, 
the Chernoff bound in (25) can be evaluated as 



E[e~ sX ] E[e sY ] 



n 



(28) 



=11 + 4s 



/>n, 



N 1 - 4s 



■N 



where 



< s < 



1 



(29) 



l+7i+Pi7i 

The tighest upper bound is obtained by selecting s that mini- 
mizes (28). This is equivalent to selecting s that maximises the 



L r 

quantity 



1 + 4s 



l+7;-P;7i 



N 



1 - 4s 



Pi7i 



l+7i+Pi7' 



■N, 



The value s that maximises the ith factor is determined by solv- 
ing the equation 



d 

rfs" 

= 



l + 4s- 



PiU 



1 + 7i - Pi% 

which gives the result 



No 



1 -4s- 



Pili 



1 + 7i + Pili 



N 



iNo' 



(30) 



(31) 



We note that this maximising value of s is independent of index 
i and falls within the allowable range of s given in (29). Using 
(31) in (28), together with (25), one finally has 



p b <E 



1 - 



Prli 
1+7* 



(32) 



The Chernoff bound in (32) can be tightened by a factor of 2 
[16, section 4.2.4], that is 



P h < 



1 L 

In 

»=i 



PHi 
1 + 7* 



(33) 



For the case of i.i.d. channels, we have p = p t = and 

7 = 7^ = 2Es j^°^ for each i = 1, 2, • ■ ■ , L. The improved 
Chernoff bound in (33) then reduces to 



Ph < 



PI 
1+7 



(34) 



which agrees with [7, eq.(4)]. Note that although the Chernoff 
bound (34) for the i.i.d. channel case can be deduced from (33) 
by setting p = pi and 7 = 7,-, the BEP given in [7, eq.(3)] can 
not be obtained from (22). This is because for the i.i.d. channels 
the denominator terms of A4 and £>; in (19) are equal to zero. 

Next, let us consider suboptimum reception. Although a nice, 
simple Chernoff bound to the error probability of the optimum 
detector can be obtained, computing the optimum bound pa- 
rameter in the case of suboptimum detection turns out to be an 
analytically cumbrous problem, especially when the order of 
diversity is large. To illustrate this point, apply Chernoff bound 
to the error probability given in (23), we have 



n'<E 



V £ |2;(fc)+2>(fc-l)| 2 



E 



«' E \Zi{k)-~Zi{k-l)\ 



n 



[1 + 4s'(l + 7l + Pili )N Q ] [l _ 4a'(l + 7i - pat)N ] 



where 



< s' < 



1 



(35) 



(36) 



4iVo(l+7i-pi7i) 

The value s' that minimizes the ith factor is determined as fol- 
lows 

[1 + 4s'(l + 7l + pai)N ] [1 - 4s'(l +7, - pi~/i)N ]} = 

(37) 

which has solution given by 



Pili 



4AM(l+ 7l ) 2 -(^) 2 ] 



(38) 



Clearly, the local value s' that minimizes each factor of (35) is 
dependent on index i, rather than equaling a constant indepen- 
dently of i. Therefore, the globally optimum bound parameter 
can only be obtained by solving 

d L 

_Jj{[l + 4 S / (l + 7 i + p i7i )iV ] x 
i=i 

[l-4s / (l+7 1 -/j i7< )JV ]}=0 (39) 



IV. NUMERICAL RESULTS 
The BEP performance is plotted in Fig.l and Fig. 2 against 
the total average received SNR per bit. The order of diver- 
sity is set to L = 2, and the fading correlation coefficient is 
set to p = 0.975 for both the identical channel case and the 
nonidentical channel case (i.e., pi = p 2 — 0.975). The ab- 
scissa represents the total average SNR per bit which is given 
by 7 b = 71 + 72- The quantity 77 is the fraction of the to- 
tal average received bit energy devoted to diversity branch 1 
(1 — T] is devoted to diversity branch 2). For example, when 
rj = 0.1, 10% of the average bit energy is devoted to branch 
1, hence, 71 = O.I7&. Clearly, in the case of i.i.d. channels, 
7 = 71 = 72 = Q.5jh- The exact BEP for the nonidentical 
channel case is plotted using (22) and (19), together with (20) 
for the optimum diversity reception and (24) for the case of sub- 
optimum diversity reception. The upper bound is plotted using 
(33) and (34) for nonidentical channel case and i.i.d. channel 
case, respectively. The exact BEP result for i.i.d. channel is 
computed using the work in [7, eq.(3a)]. 

From Fig.l and Fig. 2, several conclusions can be drawn. 
First, unequal SNR distribution among diversity branches de- 
grades the BEP performance. For instance, when the total aver- 
age SNR per bit 7b = 15 dB, the BEP is equal to 5.0234 x 10" 3 
for i.i.d. channel (7 = 71 = 72 = 12 dB), 1.065 x 10~ 2 for op- 
timum detection with r\ = 0.1 (71 = 5 dB and 72 = 14.54 dB) 
and 1.093 x 10~ 2 for suboptimum detection with 77 = 0.1. As 
■q increases to 0.5001, we can see that the BEP performance of 
optimum and suboptimum detectors for nonidentical channels 
converges to the i.i.d. channels (see Fig. 2). This observation 
also numerically validates our analytical result (22). Note that 
since the denominator terms of A; and Bi in (19) are equal to 
zero for the case of i.i.d. channels, we have set r\ — 0.5001 in- 
stead of using 77 = 0.5 for our numerical investiagtion. Second, 
for the case of nonidentical channels, in comparison with the 
suboptimum detector given in (8), the optimum detector shown 
in (6) and (7) can substantially improve the BEP performance, 
especially in the regime of high total mean SNR. For example, 
for 76 = 30 dB and r) = 0.1, the BEP is 1.616 x 10~ 3 for 
suboptimum detection and is 6.710 x 10~ 4 for optimum detec- 
tion. Moreover, as the total average SNR goes to infinity, for 
the same order of diversity L and the same value of correlation 
coefficient p, asymtotically the optimum detector for noniden- 
tical channels performs identically to the combining detector 
for i.i.d. channels. Actually, if we let 71 — > 00 and 72 —> 00 in 
(33) and 7 — > 00 in (34), it can be easily seen that the BEP floor 
of the improved upper bound (33) for nonidentical channels is 
identical to that of (34) for i.i.d. channels. 

Appendix I 

The task of the receiver is to determine from the received sig- 
nals {Si(k), Si(k — which one of two possible values 
and 7r of the phase difference A(j>(k) has maximum a pos- 
teriori probability (MAP). More precisely, the receiver will set 
A(j>(k) = nn whenever 



P 



A0(fc) 



{Si(k), Zi(k 



mti 



0,1 (40) 



which is intractable when the order of diversity is large. 



is a maximum for m 



Assuming the two data bits are equally likely, MAP detection 
is equivalent to the maximum likelihood (ML) detection. That 
is, based on {zi(k), Zi(k—l)}f =1 , we decide that A<j>(k) = irn 
whenever the PDF 



P 



{Si(k), 2i(A-l)}f =1 A0(fc) 



m = 0, 1 
(41) 

is a maximum for m = n. To simplify (41), we take the nat- 
ural logarithm for both sides and use the independent channel 
assumption, resulting in log-likelihood 

L 

log* m = ^logjp Zi(k) Zi(k- l),A(/>(k) =7rm | 
' log^p Zi(k-1) Acf)(k)=7rm | . (42) 



i=i 

L 

?:=i 



Since the third term in (42) does not affect the log-likelihood 
decision, we only need consider the second term in the compu- 
tation of the matrics. Conditioning on Zi(k — 1) and A<fr(k) = 
nm, the quantity Zj(fc) is a Gaussian random variable with 
mean 2 e^I^+ Nq h ( k - 1 ) e" m and variance 2E b R l (0)+N a - 

2E b R (o)+N ^ ]■ Using the conditional PDF of Zi(k) in (42), 
and after manipulation and simplification, one has 
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BEP COMPARISON FOR L = 2; p 



0.975. 



log* m = C + — x 

L 



Re 



{piii) 



:Zi(k) z*(k-l)e 



J7T m 



where C represents the constant term which does not affect the 
likelihood decision. The likelihood ratio test in (6) and (7) then 
follows. 
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Fig. 2 

BEP COMPARISON FOR L = 2; p = 0.975. 



